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ABSTRACT 



We report observation of the supernova remnant IC 443 (0189.1+3.0) with the 
Fermi Gamma-ray Space Telescope Large Area Telescope (LAT) in the energy band 
between 200 MeV and 50 GeV. IC 443 is a shell-type supernova remnant with mixed 
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morphology located off the outer Galactic plane where high-energy emission has been 
detected in the X-ray, GeV and TeV gamma-ray bands. Past observations suggest 
IC 443 has been interacting with surrounding interstellar matter. Proximity between 
dense shocked molecular clouds and GeV— TeV gamma-ray emission regions detected 
by EGRET, MAGIC and VERITAS suggests an interpretation that cosmic-ray (CR) 
particles are accelerated by the SNR. With the high gamma-ray statistics and broad 
energy coverage provided by the LAT, we accurately characterize the gamma-ray emis- 
sion produced by the CRs accelerated at IC 443. The emission region is extended in the 
energy band with ^gs = 0.27° ± 0.01° (stat) it 0.03° (sys) for an assumed 2-dimensional 
Gaussian profile and overlaps almost completely with the extended source region of 
VERITAS. Its centroid is displaced significantly from the known pulsar wind nebula 
(PWN) which suggests the PWN is not the major contributor in the present energy 
band. The observed spectrum changes its power-law slope continuously and continues 
smoothly to the MAGIG and VERITAS data points. The combined gamma-ray spec- 
trum (200 MeV < E < 2 TeV) is reproduced well by decays of neutral pions produced 
by a broken power-law proton spectrum with a break around 70 GeV. 

Subject headings: gamma-rays: general, supernovae: individual (IC 443) 



1. Introduction 



IC 443 is a well-studied supernova remnant (SNR), possessing strong molecular line emission 
regions that make it a case for a SNR interacting with molecular clouds. The SNR is one of the best 
candidates for revealin g the connection ani ong SNRs, molecular clouds and high-energy gamma-ray 
sources as reviewed by iTorres et al.l (j2003l ) . 



IC 443 is located in the outer Gal a,ctic plane a nd listed as a core-collapse supernova remnant 



(SNR), G189. 1-^3.0, in Green's catalog (|Greenll2004l ). The SNR has an angular extent of ~45' in the 



radio with a complex shape consisting of two h alf -shells with differ e nt radii (Shells A and B) (e.g. , 
Fesen fc KirshneiUlQSol : ISraun fc Stnm]ll986al B IPetre et al.lll988l : iRirst et al.lll99d : iLeahvl 120041 . 



Petre 



1981 



et al 



1988 



Chevaliei 



Troi? 



199S 



a et al.ll2008l) but others indicate that it i s older (20 — 30 ky) 



(e.g., 


Lozinskava 


Lee et al. 


I2OO8 





Its distance has not been measured directly but i s assumed to be ^1.5 kpc, the distance to the Gem 



OBI association to which t he SNR belongs (e.g.. lWoltierlll972l : l01bert et al.ll200ll : IWelsh k Sallmen 



2OO3I : iGaensler et al.l l2006l 'l. A pulsar wind nebula (PWN), CXOU J061705.3-h222127, has been 
found in the south ern periphery of the SNR but its assoc i ation with the SNR has not yet been 
firmlv estabhshed dKeohane et al.1ll997l : loibert et aPboOll : [Bocchino Bvkovl I2OO1I : iLeahvlbooJ : 
Gaensler et al.ll2006l : lTroja et al.ll2008l ). To this day pulsation has not been reported at the position 
of the putative pulsar. 



A general picture has been drawn from past observations and analyses th at a variety of dy- 



namical proces s es are taking pla ce in the complex structure of IC 443 (e.g., iTroia et al.l 12006 



Lee et al.ll2008l : iTroia et al.l 120081 . and references therein). The processes include: interaction of 
SNR shocks with molecular and atomic clouds of various densities which produced a break-out 
(Shell B) from Shell A as well as as sociated small-scale structur e s; interaction of the half-shells 



with another SNR G189.6+3.3 (e.g., lAsaoka Sz Aschenbachlll994l : iKeohane et al.l 119971 ): penetra- 



tion o f s hock fronts into dense molecular clouds leading to molecular line emission (e.g.. iDenover 



1979a 



Denover &: Frerkine 



Richter et al 



I995I: 



environment (Olbert et al 



Chevalier 



1981 



Huang et a l. 1986: lBurton et al. 



1988 



van Dishoeck et al. 



199 



I 



I999I : iHewitt et al.ii2006l): and interaction between the PWN and the 



2001 



Leahvll2004l : 



Gaensler et al 



2006 



Troia et al.ll2008l ). 



Of special interest for this study are the detections of high and very high energy (VHE) gamma 
rays in the IC 443 vicinity. EGR ET detected a garnma-ray source above 100 MeV, co-spatial E 



with the SNR (3EG J0617 +2238) (ISturner fc Dermer 



1995 



Albert et al 



I2OO7I 1 



Esposito et al 



1996 



Lamb fc Macomb 



I997I: iHartman et allll999l '). The MAGIC telescope discovered a VHE source, MAGIC J0616-t-225 



which is displaced with respect to the position of the EGRET source, and 
co-spatial with what appears to be the m ost massive molecu l ar cloud in the neighborhood de- 



tected in ^^CO and ^^CO emission lines ( Burton et al 



2001 



1988 



Dickman et al 



1992 



Dame et al 



Seta et all 119981 ). VERITA S has confirmed the VHE emission (VER J0616.9+2230) and re- 



solved the source to be extended (jAcciari et al.ll2009l ). The centroids of these 3 gamma-ray sources 
are displaced from that of the PWN. 

The LAT data for IC 443 provide an exciting opportunity to study the interaction of an SNR 
with the interstellar medium, cosmic-ray (CR) acceleration and subsequent injection to the Galactic 
space. The entire Milky Way has been deeply observed by the LAT and modeling of the diffuse 
emission thereon allows the emission associated with IC 443 (the "IC 443 contribution") to be 
considered separately from the underlying Galactic diffuse emission, which has contributions from 
inverse Compton scattering of CR electrons (the "Galactic IC component") and CR electron and 
proton interactions with interstellar nuclei (the "Galactic CR contribution"). In the LAT data 
the spatial extension of the contribution from IC 443 can be measured along with its broad-band 
spectrum 

This paper is organized in the following sections: A brief description of the observation, event 
reconstruction and gamma-ray selection is given in section [2j The analysis procedure is explained 
in section [3] including the instrument response function (IRF) and separation of the Galactic CR 
contribution. Galactic IC contribution, extragalactic emission and instrumental background. We 
present results on the spatial extension and spectrum of the IC 443 contribution in section [4l 
Discussion is given in section [5] and the paper is concluded in section [H 



We assume that the gamma-ray sources detected in the region are associated with locally accelerated CRs based 
on the spatial overlap with the IC 443 structure seen in the radio, IR, optical and X-ray bands. 



2. Observations, Event Reconstruction and Gamma-Ray Selection 



2.1. Observation in the Survey Mode 



The Fermi Gamma-ray Space Telescope, launched on 11 June 2008, has been surveying all 
sky with the Large Area Telescope (LAT) since August 2008. Its wide field of view (~ 2.4 sr), 
large effective area (~ 8000 cm^ at > 1 GeV), the improved point spread function (PSF) (^gg'^ ~ 
0.6° /0. 9° at E = 1 GeV) and the broad energy coverage (20 MeV — 300 GeV) bring muc h improved 
sensitivity and gamma-ray statistics over its predecessor EGRET (jAtwood et al.ll2009l ). 



The LAT is a pair-conversion telescope where a gamma-ray is converted to an e~^e~ pair. Their 
trajectories are recorded in the tracker and the subsequent shower development are sampled both 
in the tracker and calorimeter. The tracker is surrounded by a segmented anti - coincidence detecto r 
which is used to reject events induced by charged CRs (lAtwood et al.ll2009l : lAbdo et al.l l2009al ) . 
The LAT PSF 



is determined at lower energies by multiple scattering in the tracker. At 
higher energies the PSF approaches to the limit given by the granularity of the tracker channels: 
it is 0.1670.26° at 5 GeV and 0.11°/0.15° at 10 GeV (|Atwood et ahlbood v 



The LAT was operated in the nominal all-sky survey during the present observation. In the 
observation the instrument axis was tilted from the zenith toward the orbit's north and south poles 
by 35 deg or 39 deg on alternate orbits to make sky coverage uniform. The trigger rate, mostly 
on cosmic rays, was ~ 2.2 kHz in average and varied between the maximum of ~ 5.0 kHz and the 
minimum of ~ 1.6 kHz dependently on the geomagnetic cut-off rigidity. On-board filtering reduced 
the event rate to ~ 450 H z for the downlink . Data taking is disabled during passages through the 
South Atlantic Anomaly (lAbdo et al.l l2009al ) . 



2.2. Gamma-ray Selection 

Gamma-ray candidates are defined in 3 classes on the gamma-ray probability, background 
expected in orbit, current knowledge of the astronomical gamma-ray fluxes, and performance o f 



the LAT. The Diffuse class has the tightest background rejection of the three (lAtwood et al.ll2009l ). 
However the effective area becomes small and strongly dependent energy below 200 MeV. The 
averaged rate for the Diffuse class event was ~ 0.6 Hz during the observation. 

In the survey mode the Earth limb, an extremely bright source of gamma-rays, comes near the 
edge of the held of view. We have removed these gamma-rays with the reconstructed zenith angles 



^The PSF is significantly different for gamma rays detected in the front and back portions of tlie tracker described 
in tiiis section: the first and second of the two values separated by / are for those detected in the front and back, 
respectively. We use 6es and ^95 with superscripts psf, error and ext to quantify the PSF, source localization error 
and source extension, respectively. The integrated probabilities in the 2-dimensional angular radii of Oes and S95 are 
68% and 95%, respectively. For a symmetric 2D Gaussian distribution ^95 is 1.62 x Oas- 



greater than 105°. 



Cosmic-ray induced background in the Diffuse class becomes comparable in intensity to gamma- 
rays from the IC 443 region at energies below 100 — 200 MeV and above 50 — 100 GeV. The back- 
ground consists of residual cosmic rays misclassified as gamma rays and cosmic rays that convert 
in the passive material just outside of the LAT without leaving a signal in the anti-coincidence 
detector. We limit the energy range of this analysis between 200 MeV and 50 GeV where the 
effective area and the instrumental background is best undersood[§. 

The data analyzed here were obtained between 4 August 2008 and 4 July 2009. The gamma 
rays in the circular region-of-interest (ROI) of radius 15° centered at the best-fit centroid of the 
IC 443 contribution to be determined in section W7i\ {I = 189.05°, h = 3.03°) are selected for later 
ananlyses. We refer to this set of events as the data set: the key selections described here are 
summarized in Table [TJ 

Events in the data set are binned in energy at 13 logarithmic steps of 0.184 starting from 
200 MeV. The matching energy-dependent exposure is calculated based on the orbit location, point- 
ing direction, orientation and live-time accumulation of the LAT. The intensity is then calculated 
by dividing maps of counts with maps of exposure in each energy bin. 



3. Analysis Procedure 

The present analysis focuses on determination of the centroid and extension of the IC 443 con- 
tribution after separation of the Galactic CR contribution in the region. The latter will reflect the 
detailed spatial structure of the molecular clouds through pionic and bremsstrahlung interactions 
and potentially can be mistaken as a part of the IC 443 contribution. The uniform all-sky coverage 
of the LAT observation allows us to isolate the IC 443 contribution cleanly from all Galactic-scale 
contributions and determine its extension accurately. The Sourcelike analysis has been designated 
specifically for this kind of analyses. Spectral analysis has been done with the LAT Science Tool 
gtlike and has been cross-checked by Sourcelike. We give a brief description of the Instrument 
Response Function (IRF) and Sourcelike before proceeding to the fitting. 



^We are currently developing an improved event classification procedure to retain higher effective area at lower 
energies and to reduce background contaminations in the entire energy range. 

^Available from http:/ /fermi. gsfc.gov/ssc/data/access/lat/BackgroundModels. html", the Fermi Science Support 
Center url for the Science Tools. 



3.1. Instrument Response Function 



The spatial extension and spectral features of the gamma-ray emission are studied by com- 
paring the observation with predictions of source models. Predictions are made by convolving the 
spatial distribution and spectrum of the source models with the IRF and the exposure for the 
observation. 

The IRF describes the overall performance of the instrument, event reconstruction and gamma- 
ray selection. In the Fermi LAT it has b een formulated, before the launch, using an instrument 



simulation progra m (I Atwood et al.l 120091') . The simulat ion program has been calibrated against 



beam test results (jAtwood et al.ll2009l : lAbdo et al.ll2009al ) and the predicted IRF has been validated 



on several bright point sources in the early operation phase. 

The variation in the trigger rate results in variation in the fraction of the trigger-enabled time 
(the "live-time fraction") between ~ 94.3 % to ~ 81.5 %. Besides lowering the live-time fraction and 
the exposure, CR hits overlaid on a genuine gamma-ray track can reduce reconstruction efficiency 
and lead to incorrect event selection. The overall inefficiency has been found to scale linearly with 
the loss in the live-time fraction with a coefficient that depends on energy. The IRF used in the 
analysis, IRF P6_V3_Diffuse, has been corrected for inefficiency by, for example, +23 %, +16 % 
and +12 % at 200 MeV, 500 MeV and 1 GeV, respectively. 



3.2. Extension Analysis with Sourcelike 

The intensity distribution observed by LAT from the IC 443 region is shown in Fig. [T] for a 
lower (1 — 5 GeV) and a higher (5 — 50 GeV) energy bands for an area of 8° x 8° centered at {£, 
b) = (189°, 3°) with 0.1° pixelization. Spatial extension of the IC 443 contribution is determined 
on the intensity distribution using Sourcelike, an analysis tool developed by the LAT team. In 
the tool, likelihood fitting is iterated to the data set assuming spatial source models and a spatial 
background model: we use combination of a symmetric 2D Gaussian source model or a point-source 
model and the standard background model. 

The standard background model used in Sourcelike is formulated by summing the Galactic 
CR contribution. Galactic IC contribution and isotropic componenj^ given in the diffuse emission 
m odel [§: it is refer red as the background here after. All bright sources detected with the LAT 



(jAbdo et al.ll2009bl ) within 15° of the the centroid are included in the background. 



The fit is performed for the entire data in the user-determined energy range. Absolute nor- 
malization of individual background components can be constrained or unconstained in the fit: we 



^The sum of the extragalactic background, unresolved sources and instrumental background: its spatial distribu- 
tion is assumed to be isotropic. 

^ glLiem_v02.fit and isotropic_iem_v02.txt available from the url given in footnote 4. 



leave the diffuse emission model as one unconstained component and so are all bright sources in the 
ROI in the Sourcelike fit. The difference in Test Statistic (TS) values between the best-fit Gaussian 
distribution and the best-fit point-source which is 2A/o5(Likelihood) gives a measure of statistical 
significance of the extension. We refer to this difference as TSext in this paper. 

The ROI is energy dependent in the Sourcelike: 15° at 200 MeV and shrinks to a minimum 
of 3.5° at 50 GeV, which is at least a factor of 20 larger than ^gg'^ of the LAT at the same energy 
and more than a factor of 10 larger than the spatial extension (^gf*) of the source to be determined 
later. 

4. The IC 443 contribution 

4.1. Spatial Extension of the IC 443 Contribution 

Two energy bands, 1 GeV< E <5 GeV (the lower energy band) and 5 GeV< E <50 GeV (the 
higher energy band) have been selected to study the spatial extension of the IC 443 contribution. 
Sourcelike has been run for events in the two energy bands separately as well as in the combined 
energy band under a 2D Gaussian and point source hypotheses. The best-fit results are summarized 
in Table [2l 

The number of gamma rays in the fitted Gaussian distribution is 4972 for 200 MeV< E <1 GeV, 
1597 for 1 GeV< E <5 GeV and 236 for 5 GeV< E <50 GeV. For a given PSF, the accuracy of cen- 
troid determination is predicted to improve proportionally to the inverse of the square-root of the 
number of events. The accuracy quoted in Table [2] is consistent with this prediction for the effective 
PSF averaged over events in the energy bands 1 GeV< E <5 GeV and 5 GeV< E <50 GeV. 

The difference in TS (TSgxt) between the symmetric 2D Gaussian and point hypotheses is 
+106 to +121 (10.3 to 11.0 0-0 for the 1 - 5 GeV band and -^212 to -^81 (14.6 to 9.0a) for 
the 5 — 50 GeV band. The centroids for the two bands are consistent within 0.04°. The IG 443 
contribution is extended to 0g|* = 0.26° — 0.27° in the two energy bands. The centroid in the high 
energy band is displaced southwards by ~ 0.04° (~ 1.5cr) in the Galactic coordinate from that in 
the low energy band. 

The results on source location and extension are robust: TS values have been examined at 
discrete points offset from the best-fit location and extension to confirm the fit. To verify the fit 
further, we have generated 100 simulated sets of events assuming the best-fit centroid, extension and 
background with the LAT Science Tool gtobssim. The simulated data are then processed through 
Sourcelike under 2D Gaussian and point source hypotheses. The distribution of TSgxt between the 
two hypotheses is consistent with the values given in Table O 



^The two TSea:t valuGS quoted are: the first one for that used in P6_V3_difi'use; and the second one for a worst-case 
PSF to be used later in section [S] to obtain a conservative systematic error. 



The point source hypothesis is rejected at TSea;t> 81 or > 9 o" independently in the two energy 
bands and at TSea;t> 212 or > 14 a in the combined energy band. The extensions in the two energy 
bands are mutually consistent within the errors given in Table [2j 

The radial profiles of event distribution around the centroid is shown in Fig. [2] for the low and 
high energy bands together with the profile predicted for the point source hypothesis which is the 
LAT PSF weighted with the spectral distribution of the events Sourcelike has associated with the 
source under the point-source hypothesis. 

Extension was poorl y determined for E= 200 MeV — 1 GeV because of the large PSF of the 



LAT in the energy range (jAtwood et al.ll2009l ). However, the centroid and extension are consistent 
with the extension determined above 1 GeV and given in Table [21 Hence we assume the same 2D 
Gaussian distibution in the entire energy range. 



4.2. Spectrum of the IC 443 Contribution 



The spectrum of the IC 443 contribution is fitted by the Science Tool gtlike, the Fermi standard 
tool, as well as by Sourcelike. In gtlike, we have to assume a spatial template for all spectral 
components included in the fitting. The data set is assumed to be a sum of three contributions: 
the best-fit 2D Gaussian distribution given for the E = 1 — 50 GeV range in Table[2]which represents 
the IC 443 contribution the background whose spatial distribution is represented by the sum of 
glLiem_v02.fit and isotropic-iem-v02.txt; and the bright sources listed in lAbdo et al. ( 2009b ) in the 
square region of 8° x 8° centered at the best-fit centroid {£ = 189.05, b = 3.03). 

The fitted IC 443 spectra from gtlike and Sourcelike agree well within the total error. We 
adopt the spectrum obtained with gtlike and tabulate in Table [3l It is converted to the spectral 
energy density (SED) and shown by circles with error bars in Fig. [3l Squares with error bars in the 
figure is the background spectrum normalized to the solid angle subtended by 6*9!* = 0.45° around 
the centroid given in Table [2j The IC 443 contribution is approximately 20 times higher than the 
background in the entire energy band. 

The SED of the IC 443 contribution thus determined has been fitted with a single power-law 
and broken-power-law models: the results are tabulated in Tabled The single power-law fit fails 
to represent the spectrum giving a large reduced chi-square of ~ 9 while the broken power-law 
fit represents the overall shape quite well giving a small reduced chi-square (~ 1.0) as shown in 
Table [H The SED of the IC 443 contri bution is plotted w ith those from previous observatio ns, 
EGRET (|Esposito et allll99d v MAGIC (lAlbert et al.ll2007l ;i and VEB I TAS ^ Acciini et ahlbood l in 
Fig.H 



The extension could not be determined at a high statistical significance for E = 200 MeV —1 GeV but the spatial 
distribution of gamma rays is consistent with those given in Tabled 



4.3. Systematic error in determination of the centroid, extension and spectrum 



When we determine the centroid of the IC443 contribution, uncertainty in the spatial distri- 
bution of the Galactic diffuse emission adds to the systematic error. The spatial template is taken 
from the standard diffuse emission model, glLiem_v02.fit. To confirm our analysis, we have fitted 
the data set with the standard version of GALPRO P for Fermi LAT (GALDEF 54 59Xvarh8S) 



Strong fc Moskalenkdll998l: IStrong et a 



200 ll) and H I survey b y 



Kalberla et al 



2001 , 12OO9I ) made of the CO line survey by I Dame et al 



(I2OO5I') as well as with a gas model made of the ™ap 



by iDobashi et al.l (|2005l ) and H I survey by iKalberla et al.l (j2005l ) . The two alternate gas models 
have given centroids consistent with that given in Table [2l 

The residual misalignment of the LAT and the star tracker can also contribute to the systematic 
error: the source localization has been verified on orbit using bright point sources to ±30 arc-sec as 
of August 2009. The combination of all errors described here gives the overall systematic localization 
error in Table O 

Our flux measurement depends on the knowledge on the effective function of gamma- 

ray energy. We estimate systematic error in the effective area to be 10 %, 5 %, and 20 % at 
E=100 MeV, 562 MeV and > 10 GeV respectively 

Uncertainty in the background used in gtlike and Sourcelike can also introduce error in the flux 
measurement. This uncert ainty is estimated to be ~ 20% for 200 MeV — 1 GeV and ~ 30% for > 
1 GeV of the background dAbdo et al.ll2009cl ldlv Systematic error at each energy bin is determined 
through a linear interpolation in logio{E) among the values quoted above. 

The PSF used in this analysis has been de rived on the detector simulation which was itself 
verified in accelerator tests (jAtwood et al.ll2009l ). As gamma-ray statistics improves, the PSF will 
be updated against measurement on bright point sources. In the present study, we have used a 
preliminary upper limit to assess possible systematic error introduced by inaccurate formulation 
of PSF: we consider this as the "worst-case" PSF. The worst-case PSF (^gg'^) gives a widest limit 
while the for E > 5 GeV is about 40% larger than the default PSF. The source centroid comes out 
to be consistent within the total error when Sourcelike is run with the worst-case PSF. We have 
included the difference in the systematic error given in Tables [2] and [5l 



5. Discussion 



The IC 443 system consis ts of a complex d i stribution of ru olecular and atomic clouds in the 
southern rim of Shell A (e.g., ISnell et al.l l2005l : iLee et al.l l2008l . and references therein). Molec- 
ular clouds wrap around the southern rim and the bound a ry region between Shells A and B 



(Beichman et al 



19881 : ISeta et al.1 Il998l : lOickman et al.1 Il992l : lOame et al.l I2OOII ) . Molecular lines 



from shocked gas have been found in several clouds suggesting interaction with the blast-wave at 



multiple sites (e.g.. ICornett et al.lll977l : iDenover &: Frerking]|l98ll : iHuang et al.lll986l : iBurton et al 



19881 : bickman et~allll992l : Ivan Dishoeck et alJEggj : iRichter et al]|l995l : ISeta et alll998l : ISnell et al 



2005 ) . A prominent band of HI ga,s has also been found in the southeastern part of Shell A (|Denovei 
19781 : iGiovanelli &: Havneslll979l: iLee et al.ll2008l. and references the rein.) . Some parts of the H I 
gas are found to be shocked (jBraun Stromlll986al : iLee et al.ll2008l ). 



The total mass of the molecular gas in the re gion is estirnated to be ~ 1 x 10^ Mq ( Torres et al 



20031 ). of which only a small fraction is shocked (jSnell et al.ll2005l). The t otal mass in the H I belt 
is estimated at ~ 730 Mq of which ~ 500 Mq is shocked (ILee et al.ll2008l ). Despite past extensive 
observations and analyses, little is known about how the multiple shell-like structures are spatially 
correlated and where one or more supernova explosions took place. 

In the group of shocked molecular clouds schema tically showii in Fig. El Cloud G@ lies closest to 



the centroid of the MAGIC and VERITAS sources (IHuang et al, 
by_~ 8^ and overlaps with a non- shocked C O cloud formation 



198d). It appears to be extended 



Huang et al.l Il986l : iBurton et al 



19881 : iDickman et al.lll992l ). IChevalieij (j 19991 ) has suggested that interaction between Shell A and 
Cloud G is responsible for the gamma-ray emission observed by EGRET. Existence of an OH 
Maser in t he cloud suggests that the densit ies reach ^ 10^ cm ~^ ( Frail et al.lll996l : iHewitt et al 



20061 . 12OO8I ). A later CO line observation bv ISnell et al.l (j2005l ) found a compact core of extension 
~ 1' — 2' in Cloud G at (i = 189.03°, b = 2.90°). Fainter maser emission has also been found in 
Clouds B a nd D at = 189.18°, b = 2.97°) and {£ = 189.25°, b = 3.13°), respectively (|Hewitt et al 



2006 



20081 ). 



The locations an d extensio ns of the gainma-ra y emission fr om IC 443 detected by EGRET 
(|Hartman et all 1 19991 ). MA GIG dAlbert et al.ll2007l ). VERITAS (|Acciari et aklbood ). and Fermi 
LAT are summarized in Tableland shown in Fig. [5j Fermi LAT gives the best source localization 
or the smallest error circle (^gg'^"'^) for the 2D Gaussian centroid (0.03°) and a precise determination 
of the source extension (0.27° it 0.01° (stat) it 0.03° (sys)). Our centroid determined for 1 GeV < 
< 50 GeV, {i = 189.05°, b = 3.03°), is 0.05 deg away from the EGRET source (3EG J0617+2238) 
but within 0'^°^ of their localization error; 0.15 deg from the MAGIG source (J0610+225) which 
is at more than 5 times their localization error {0'^°^); and 0.12 deg from the VERITAS source 
(VER J0616. 9+2230) or at 1.5 times their localization error {6\ 



errors 
68 J 



The measured source extension, = 0.27° =b 0.01° (stat) it 0.03° (sys), is comparable with 



eif = 0.24° ± 0.05° (stat) ± 0.06° {sysj^ given by VERITAS. The two extended regions overlap 
almost completely. The three shocked clouds with OH maser (Clouds B, D and G) are within our 
measured extension and so are ot her shocked c l ouds (Clou ds C, E, F and H). T he PWN localized 
at {£ = 189.227°, b = 2.897°) bv loibert etlo] (|200ll ) and ICaensler et all koo4 ) is 0.26 deg away 



from our centroid but within our measured extension ^gf*. 



The Fermi spectrum of the IC 443 contribution shown in Fig. [3] is flat between a few 100 MeV 



'Labeling is given bv lPenovei] l|l979bl ): iHuang et all (|l986h . 

We assume the extension is modeled by a symmetric 2D Gaussian and converted to ^$8 ~ 1-51 x 9ia-. 



and ~ 3 GeV, suggesting the origin being mostly neutral pions produced by protons 0. The dashed 
line in the figure represents the gamma-ray spectrum expected from a 10^ M.q cloud bombarded 
with the Galactic CR protons predicted at IC 443 scaled up by a f actor of 100. The Ga.lactic 



CR spectrurn is taken fr om the standard GALPROP (54_59Xvarh8S) (IStrong &: Moskalenkq 



1998 



Strong et al.l l200ll . |2009|) and the parameterized cross-section for ^ 7 by iKamae et al.l (|200Q ) . 
In GALPROP the Galactic CR proton spectrum depends on the radius from the Galactic Center 
and the displacement from the Galactic Plane. The spectrum at the radius of IC 443 is ~ 10 % 
lower than that in the solar vicinity and has a power-law shape with index ~ 2.7. 

The Fermi SED is compared with those of EGRET MAGIC and VERITAS in Fig. H The 
EGRET spectrum is consistent with our spectrum except for their 3 GeV point. MAGIG and 
VERITAS do not overlap with the LAT in the energy coverage. Their fluxes are consistent with 
ours if extrapolated down to ~ 50 GeV by assuming their measured power-law indices. 

Since the source regions of Fermi LAT and VERITAS overlap within their respective uncer- 
tainties listed in Tabled we can judiciously proceed to fit the 2 spectra with one spectral model. 

On the assumption that the distance is d = 1.5 kpc, the isotropic luminosity of IC 443 inte- 
grated over the energy band (0.2 — 50 GeV) is 1.2 x 10^^ erg/s. Electron bremsstrahlung can hardly 
explain the observed IC 443 gamma-ray emissivity: the cross-sections for bremsstrahlung and pionic 
gamma-ray emission are similar in the present energy band, so the bremsstrahlung-to-pion ratio is 
approximately the ratio of CR electron and proton fluxes which is ~ 0.01. The observed gamma-ray 
flux is too high for bremsstrahlung to be the dominant process. Inverse Compton scattering can 
not explain the observed IC 443 gamma-ray emission either: the gas density of the emission region 
is ~ 50 — 100 cm~^ and the Compton-to-bremsstrahlung ratio is ~ 0.01 — 0.001 for the seed photon 
density of the cosmic microwave background. We note that there is no bright source of seed photons 
known in the region of the IC 443 contribution. The gamma-ray energy will be strongly bound by 
the electron spectrum which likely rolls down similarly as the proton spectrum. Our observation, 
however, does not rule out a small contribution from bremsstrahlung near the minimum of the 
present energy band. 

In a hadronic scenario, the observed photon spectrum up to TeV energies can be well fitted 
by an underlying pion-producing proton population with a broken power-law spectrum Fp(Tp) = 
5.9 X 10-2 (Tp/69 GeV)-° (lO'^Me/Mga^) cm-2s-iGeV-\ where a is 2.09 ±0.04 for Tp < 69 GeV 
and 2.87ib0.07 for Tp > 69 GeV respectively, and Mgas is the gas mass in the interaction region. The 
error (statistical) in the fitted break energy is ±25 GeV and the chi-square for the best-fit broken 
power-law model is 9.6/14 per degree-of-freedom. Assuming the gas density (n cm"^) is uniform and 



'We include cosmic-rays and target nuclei heavier than the proton (the alpha particle and heavier nuclei) in protons 



throughout this paper. In the approximation we adopt here (|Gaisser fc Schaefer lliggd ). these can be accounted for 
by multiplying a nuclear factor (~ 1.7) without changing the CR proton spectrum. The known gamma-ray producing 



particle processes which do not go through neutral pions (e.g.. 
less than 1 % in the present energy range. 



77 and direct photon production) contribute 



the proton spectrum is the broken power-law everywhere in the interaction region, the total energy 
of the interacting protons is given by Wp(> 0.5 GeV) = 5.6 x 10^^(n/240 cm~^)~^(d/1.5 kpc)^ erg. 



gas 



Note that the pion production threshold is ~ 0.5 GeV. Taking M, 
ranging from 47r/3x (0gg*d)^ = 5 x 10^^ cm^ up to 47r/3x (^gf^d)^ 



10 M0 and a gas volume 
2 X 10^^ cm^, we obtain n 



~ 60 — 240 cm~^, giving Wp(> 0.5 GeV) = (0.56 — 2.2) x 10^^ erg. We note that energies carried 
by local nuclear cosmic rays outside of the interaction region and by local leptonic cosmic rays are 
not included in the estimation. The fitted gamma-ray spectrum is shown in Fig. [3] and Fig. [H We 
note that inclusion of the MAGIC points in the fit does not change the above results. 

Broad-band gamma-ray spectral models ha ve be en proposed as s uming CR interaction with 
interstellar gas in IC 443 by iTorres et al.1 (|2008l ) and Izhang Fa^ (|2008l ). iTorres et al.1 (|2008l ) 
model the CR diffusion in the SNR to allow spectral differenc es in the protons inte racting with 
the ambient gas in the shell and in a detached molecular cloud. IZhang fc Fang! (j2008l ) predict one 
contribution to come from the SNR shell evolving in the interstellar matter and the other from CR 
i ntera ction with molecular clouds. One difference between the two models is that IZhang &: Fang 
(|2008l ) include inverse Compton scattering as a possible emission mechanism. 



The combined spectrum of Fermi and VERITAS gives a strong constraint to spectral models 
for the IC 443 contribution. Since the spectrum of the dominant CR component (proton) is rolling 
over at ~ 70 GeV, secondary electrons and positrons can only contri b ute a t energies < 7 GeV. 
This constrains the parameter space of the model by Zhang fc Fang (j2008l ). The overlap between 
the Fermi and VERITAS spatial extensions and the smo oth spectral transi tion from Fermi to 
VERI TAS constrain the parameter spaces of the models by lTorres et al.l (120081 ) and lZhang fc: Fang 

hood ). 



We discuss briefly about possible mechanisms behind the broken power-law form of the proton 
spectrum deduced from the gamma-ray observation of IC 443. The most obvious one is escape of 
highest energy CRs from the acceleration site. When accelerated CR protons exceed the maximum 
energy determined by the magnetic field and linear-scale of the aceleration site, they escape into 
the Galactic space. Theory of diffusive shock acceleration (DSA) assumes spherically symmetric 
morphology and predicts the CR spectrum to roll over exponentiall y at the maxmum energy. The 
maximum energy depends on the condition of the acceleration site: iPtuskin &: Zirakashvilil (j2005l ) 
have incorporated various instabilities in DSA and predict the maximum energy as a function of the 
SNR age. For the age of IC 443 (^ 30 kyr), the maximum proton energy can be around 100 GeV 
(Fig.l of iPtuskin &: Zirakashvilil l2005l ) or near the observed break energy (~ 69 GeV) beyond 
which the spectrum is assumed to cut-off exponentially. We have fitted the observed gamma-ray 
spectrum with a single power-law, exponentially cut off proton spectrum to get chi-square per 
degree-of-freedom of 30.3/15 much higher than 9.6/14 for the broken power-law spectrum. This 
simple statistical test therefore suggests that the observed broad-band gamma-ray spectrum is 
inconsistent with the simple DSA-based CR escape scenario which predicts a simple powerlaw with 
an exponential cutoff in the proton spectrum. We also note that DSA has mostly been studied for 
uniform gas densities around lcm-3, while the gas around IC443 is inhomogenous and in various 



shocked molecular clouds around the remnant also denser. 



Historically, after the discoverie s of the EGRET sour ce (jSturner &: Dermeiill995l : lEsposito et al 



lOQa) and the har d X-ray source (IKeohane et al.l 119971) b ut before the discoveries of the PWN 
(jOlbert et al.ll200ll ) and MAGIC source (lAlbert et al.ll2007l ). models have been propos ed to explain 



the e mission between ~ 5 keV and ~ 5 GeV by bremsstrahluiig in d ense clouds (e.g., iBykov et al 



2000) with possible mix of synchrotron (e.g., ISturner et al.lll997l ). The SNR wa.s also studied 
as a part of non- linear shock evolution in various environments (e.g., iBaring et al.lll999l ). While 
their predictions for IC 443 are not supported by the later observations including the present one, 
parameters in these studies can be readjusted to describe the bremsstrahlung contribution discussed 
below. 

The bremsstrahlung likely makes a non-negligible contribution below = 200 MeV where 
the EGRET data points exceed the best-fit pionic spectrum (see Fig. As our understanding of 
the IRF and cosmic-ray-induced background improve, analysis will be extended to energies lower 
than 200 MeV and the bremsstrahlung spe ctrum component will be determined accurately. The 



hard X-ray SE P measured by Beppo-S AX (jBocchino Bvkovl l200d ) is substantially higher than 
that by XMM ( Bocchino &: Bvkovll200ll ) : which may s u ggest bre msstrahlung contribu tion near the 



PWN location as has been discussed by ISturner et al.l (jl997l ) and lBykov et all ()2000l ). 



6. Conclusions 

We have studied gamma-ray emission from the nearby SNR IC 443 (G189. 1+3.0) using the 
first 11 months of science data from Fermi LAT. The uniform sky coverage and high gamma-ray 
statistics of the observation have enabled us to separate the genuine IC 443 contribution from the 
emissions by Galactic CRs on interstellar gas, inverse Compton scattering by Galactic CR electrons 
on large-scale interstellar radiation field, extragalactic sources and instrumental background. 

Based on the extension study described in subsections 14.11 and the spectral analysis described 
in subsection 14.21 as well as discussions given in section [Sj we conclude that: 

• The gamma-ray emission from IC 443 is detected at ~ 860" level: the emission is extended with 
68% containment angular radius 0g|* = 0.27°±0.01°ib0.03° in the energy range between 1 GeV 
and 50 GeV. The extension remains unchanged within error in the low (1 GeV < E < 5 GeV) 
and high (5 GeV < E < 50 GeV) energy bands. 

• The centroid of the emission moves at ~ 1 — 1.5(7 level toward that of the VERITAS source as 
the energy band changes from 1 — 5 GeV to 5 — 50 GeV. The centroid is inconsistent with the 
PWN location, suggesting that the PWN is not the major contributor in the present energy 
range. 

• The centroid of the emission is consistent with that of EGRET (3EG J0617+2238), dis- 



placed more than 5 x 9^1'''^ (MAGIC) from that of MAGIC (J0610+225), and at 1.5 x 
e^l'''^ {VERITAS) that of VERITAS (VER J0616.9+2230). 

• The extended source region overlaps almost completely that of VERITAS. A group of molecu- 
lar clouds (Clouds B, C, D, F, and G), the SNR shell and the PWN arc within the overlapping 
region (^gg*), leaving possibility that some or all of them contribute to the observed emission. 

• The SED can not be represented by a single power-law but is consistent with a broken power- 
law with a break at = 3.25 ± 0.6 GeV. 

• The SED has a broad peak between a few 100 MeV and ~ 5 GeV which is consistent with the 
majority of the emission coming from neutral pion decays. For the emission being hadronic 
originating from a single proton population, the underlying proton spectrum is consistent 
with a broken power-law shape (chi-square per degree-of-freedom of 9.6/14) but not with an 
exponential cut-off (30.3/15). For the estimated total mass of interacting gas of 10^ M©, the 
total energy in the pion-producing protons is estimated to be (0.56 — 2.2) x 10^^ erg. 

Higher statistics is needed to establish association or non-association of the gamma-ray emis- 
sion with the molccTilar clouds and/or the PWN as well as CR injection process from the SNR into 
the Galactic space. Identification of the emission mechanisms and underlying CR spectra effective 
in individual sites will follow after such studies. 

Fermi LAT is expected to accumulate needed statistics well within the planned mission lifetime. 
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Table 1. Selection for the Data Set 



Parameter 

Time Period (MET) 

Energy Range 

ROI 

Photon Class 

Additional Cut 



Value 



239557417 - 268416079 
200 MeV - 50 GeV 
< 15° in radius 
Pass 6 Diffuse 
Zenith angle < 105° 



Table 2. Centroid and Extension of the IC 443 Contribution 



Model i (°) h (°) e^^^°^ (°) ^ e^' (°) ^ eif (°) ^ TS^^t 

1 - 5 GeV 

Point Source 189.05 3.04 0.02 ... ... 

Gaussian 189.05 3.05 0.02 0.27 ± 0.03 0.44 ± 0.04 +106/ + 121 

5-50 GeV 

Point Source 189.05 2.98 0.03 ... ... 

Gaussian 189.06 3.00 0.03 0.26 ± 0.04 0.42 ±0.07 +81/ + 212 
1 - 50 GeV 

Point Source 189.05 3.02 0.02 ... ... 

Gaussian 189.05 3.03 0.02 0.27 ± 0.03 0.45 ± 0.05 +212/ + 362 



'^Errors of centroids and extensions quoted include systematic errors. 

^The two values shown in this column are the lower and higher of the TSext obtained with 
the default PSF (the first) and worst-case PSF (the second) described in the text. 



Table 3. Spectrum of the IC 443 Contribution 



-^center [JVleVJ 


aiy/aiii [cm s MeV J 


aiy lahi (stat errorj 


dN/dE (sys error) 


247.31 


8.63xl0~^° 


4.80x10-11 


6.37x10-11 


378.19 


4.11x10"^° 


1.87x10-11 


2.52x10-11 


578.32 


1.72x10"^° 


7.81x10-12 


8.89x10-12 


884.36 


8.31x10-^1 


3.59x10-12 


6.11x10-12 


1352.34 


3.25x10-11 


1.65x10-12 


3.12x10-12 


2067.96 


1.51x10-11 


8.36x10-13 


1.78x10-12 


3162.27 


5.81x10-12 


4.06x10-13 


8.14x10-13 


4835.67 


2.40x10-12 


2.05x10-13 


3.90x10-13 


7394.58 


7.49x10-13 


8.98x10-1^ 


1.38x10-13 


11307.60 


2.80x10-13 


4.29x10-1^ 


5.61x10-1^ 


17291.30 


9.13x10-1^ 


1.93x10-1^ 


1.82x10-1^ 


26441.41 


2.95x10-1^ 


8.94x10-15 


5.90x10-15 


40433.51 


6.46x10-15 


3.35x10-15 


1.29x10-15 



Table 4. Spectral Fit to the IC 443 Contribution 





Model 


71 


72 


Efcreofe (GeV) 


F2oo^ (10-^ cm-2 s-i) 


xVdof 


IC 443 


Broken PL 


1.93 ± 0.03 


2.56 ± 0.11 


3.25 ± 0.6 


2.85 ±0.07 


8.9/9 


IC 443 


PL 


2.08 ± 0.02 






3.00 ±0.07 


90/11 



^Total flux integrated above 200 MeV obtained with gtlike assuming the best-fit broken power-law 
model and the best-fit 2D Gaussian spatial distribution. 



Tabic 5. STimmary of Locations and Extensions of the Gamma-Ray SoTirccs 



Observation 


in 


bn 


of localization (°) 


^gg* of extension (°) 


EGRET 


189.00 


3.05 


0.13 {911'^"'') 


N/A 


MAGIC 


189.03 


2.90 


± 0.025 (stat) ± 0.017 (sys) 


N/A 


VERITAS 


189.07 


2.92 


± 0.03 (stat) ± 0.08 (sys) 


0.24 ± 0.05 (stat) ±0.06 (sys) 


Fermi 


189.05 


3.03 


± 0.01 (stat) ± 0.02 (sys) 


0.27 ± 0.01 (stat) ± 0.03 (sys) 



Lldeg] 




Fig. 1. — Intensity map of the IC 443 region in the 1 — 5 GeV (left) and 5 — 50 GeV (right) bands. 
Units of intensity are 10~^ cm~ ^ s~^ sr~^ f o r the color scale. The overlay is the 2.7 GHz radio 
continuum contours taken from iFurst et al.l (|l99Cl ). The insets are the spectrum- weighted LAT 



PSF for each energy band, with the white circles showing the corresponding 0\ 



68 • 




Fig. 2. — Radial profile of the gamma rays SourceLike associated with the IC 443 contribution 
(data points), compared with the effective spectrum- weighted LAT PSF (solid red line) and the 
fitted 2D Gaussian extended model (dashed blue line). The left and right panels correspond to 
the low (1 — 5 GeV) and high (5 — 50 GeV) energy bands, respectively. The points in the count 
profile are plotted at the weighted average radial positions within their respective bins. The vertical 
dash-dot and dashed lines correspond to the fitted extension 0q|* and given in Table [2] for the 
respective energy bands. 
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Fig. 3. — The gamma-ray spectrum of the IC 443 contribution: The upper and lower data points 
represent the IC 443 contribution and the total background, respectively. The background has been 
scaled to match the solid angle subtended by a disk of radius 0g|* = 0.45 deg. Errors are shown 
by the bars (statistical) and the grey band (systematic) . The lines represent the pionic gamma-ray 
spectra produced by the Galactic CR proton at IC 443 scaled up by a factor of 100 (dashed) and 
by the locally accelerated proton population with the best-fit broken power-law spectrum described 
in the text (dot-dash). 
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Fig. 4. — Broadband spectral energy density of the 4 sources: EGRET (purple triangles), MAGIC 
(blue squares), VERITAS (red diamonds) and Fermi (black circles). The solid line is the same as 
the dot-dash line in FiglS) The systematic and statistical errors of the Fermi data points are also 
the same as in Fig. [3l 
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Fig. 5. — Locations and extensions of the 4 gamma-ray sources: EGRET centroid (A); MAGIC 
centroid (v)j VERITAS centroid (star) and Fermi LAT centroid (o). The respective localization 
errors as tabulated in Table [5] are shown as crosses. Best-fit spatial extensions of the Fermi (cross- 
hatched band) and VERITAS (striped green band) sources are drawn as rings with radii of ^gf* 
and widths of itlcr error. The PWN location is shown as a dot. Contours are the locations and 



shapes of the local shocked molecular clouds taken from iHuang et al.l (|l986l ). 



Fig. 1. — Intensity map of the IC 443 region in the 1 — 5 GeV (left) and 5 — 50 GeV (right) bands. 
Units of intensity are 10~^ cm~ ^ s~^ sr~^ f o r the color scale. The overlay is the 2.7 GHz radio 
continuum contours taken from iFurst et al.l (|l99Cl ). The insets are the spectrum- weighted LAT 



PSF for each energy band, with the white circles showing the corresponding 0\ 
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